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ABSTRACT 

Context. Recent near-infrared polarization measurements of Sgr A* show that its emission is significantly polarized during flares and consists 
of a non- or weakly polarized main flare with highly polarized sub-flares. The flare activity suggests a quasi-periodicity of ~20 minutes in 
agreement with previous observations. 

Aims. By simultaneous fitting of the lightcurve fluctuations and the time- variable polarization angle, we address the question of whether these 
changes are consistent with a simple hot spot/ring model, in which the interplay of relativistic effects plays the major role, or whether some 
more complex dependency of the intrinsic emissivity is required. 

Methods. We discuss the significance of the 20min peak in the periodogram of a flare from 2003. We consider all general relativistic effects 
that imprint on the polarization degree and angle and fit the recent polarimetric data, assuming that the synchrotron mechanism is responsible 
for the intrinsic polarization and considering two different magnetic field configurations. 

Results. Within the quality of the available data, we think that the model of a single spot in addition to an underlying ring is favoured. In this 
model the broad near-infrared flares of Sgr A* are due to a sound wave that travels around the MBH once while the sub-flares, superimposed 
on the broad flare, are due to transiently heated and accelerated electrons which can be modeled as a plasma blob. Within this model it turns 
out that a strong statement about the spin parameter is difficult to achieve, while the inclination can be constrained to values i 35° on a 3<x level. 
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1. Introduction 

It is by now widely accepted that the compact radio source 
Sagittarius A* (Sgr A*) is the manifestation of a massive (M 
= 3.6 ■ 1O 6 M ) black hole (MBH) at the center of the Milky 
Way (e.g. Sc hodel et al. EDDll Genzel et aL EOOH Ghez 
et al. l2005al Eckart et al. l2006al b'). It is the closest MBH, with 
a distance of only ~ 7.6 kpc (Eisenhauer et al. 2005 1. The first 
observation of its highly variable near-infrared (NIR) coun- 
terpart has been reported by Genzel et al. (2003). The short 
rise and fall timescales of these NIR flares point to physical 
processes within a few Schwarzschild radii (Rs) of the MBH. 
Simultaneous NIR/X-Ray observations showed that the high 
variability at these wavelengths have a common physical ori- 
gin, with a spectrum pointing at synchrotron and synchrotron 
self-Compton models with rapid cooling of transiently heated 
electrons (e.g. Eckart et al. I2004I l2006al Liu et al.[2006 Liu 
& Melia 2002 Yuan et aL EOoHEOOH Markoff et al. EOOlV 

In a few NIR and X-Ray flares a suggestive quasi- 
periodicity of roughly ~ 20 minutes has been reported (Genzel 



et al. EOoH Eckart et al.l2006b Belanger et al. 2006 ). It man- 
ifests itself in the lightcurve as sub-flares superimposed on the 
underlying main flare (see Figure 0. It is not clear yet, how 
these quasi-periodicities (QPOs) are created. The timescale of 
the lightcurve variations and the rather small volume around 
Sgr A*, where these variations originate, indicate that rapid 
motion in strong gravity is involved. Frequencies in micro- 
quasar QPOs scale with the inverse mass of the BH, and the 
frequency of the QPOs in Sgr A* indicates that an extrapo- 
lation from the stellar-mass BH case in microquasars to the 
MBH case is applicable (Abramowicz et al. 2004 McClintock 
& Remillard 2004 and references therein). 

The timescale of ~ 20 minutes is comparable to the orbital 
timescale T near the innermost stable circular orbit (ISCO) of 
a spinning MBH (Bardeen et al. ll972l . This can be written as 



T = 110 (r 3/2 + a*) 



M 



3.6 • 10 6 M Q 



[sec], 



(1) 



where a* is the BH dimensionless spin parameter (-1 < 
a* < 1), and r is a circumferential radius within the equato- 
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rial plane given in units of the gravitational radius r„ = GM/c 2 . 
Therefore it appears reasonable that corresponding frequencies 
are present in the flare lightcurves. In this context, we would 
also like to point out the rapid variability in the light curve of 
the NIR flare of 16 June 2003 shown in Figure 1 (originally 
reported by Genzel et al. l2003> . Variations greater than 5<x can 
be seen on timescales of a few minutes, inferring an upper limit 
of the size of the source of less than 10 Schwarzschild radii if 
assuming that the cause responsible for the flare propagates at 
the speed of light. It can be expected, however, that the actual 
signal in the source propagates with a speed comparable to the 
sound speed or to the Alfven speed, i.e. a few orders of mag- 
nitude slower than c. Therefore it appears again reasonable to 
assume orbital motion of a compact source as the cause of the 
rapid variability. 

Furthermore, global MHD simulations show that the in- 
ner accretion flow is very inhomogenous and nonuniform, so 
that local overdensities build up (Machida & Matsumoto 2003 
DeVilliers et al. 2003 1. Modeling these overdensities as hot 
spots orbiting around the BH give lightcurves very similar to 
the observed ones (Broderick & Loeb 2005, 2006; Eckart et 
al. 2006b I. Recent polarimetric measurements showed that the 
sub-flares have a degree of linear polarization of up to ~ 20% 
and that they come along with a swing in the polarization an- 
gle (Eckart et al. 2006b I. The new information provided by the 
polarization state is extremely useful because it breaks the de- 
generacy of various model parameters. In this paper we demon- 
strate that the available data are consistent with an orbiting spot 
model and that they offer a new avenue to constrain BH param- 
eters in the future. 

In the next section, we first re-analyse the NIR flare 
observed by Genzel et al. J2003t . Using the method of 
Vaughan (2005 1 we examine the significance of the peak in the 
periodogram at ~20 min. After that, we address the question 
what the nature of the broad underlying main flare might be. 
We then use this, together with the hot spot model for the sub- 
flares, to fit the recent NIR polarimetric observations by Eckart 
et al. (2006b I. Details of the model and the fitting procedure 
are given. In Sect. |4] we draw conclusions and discuss caveats, 
with special emphasis on the effects and analytic properties of 
higher order images. 

2. Quasi-periodicity in NIR flares from Sgr A*? 

Up to now, quasi-periodicities in NIR flares from Sgr A* have 
been reported just for the three K-band flares described in 
Genzel et al. |2003 ) and Eckart et al. d2006bl l. observed at the 
VLT. Their detection at longer wavelengths has not yet been 
reported on. The reason for this may be indicative for a lower 
sub-flare contrast at longer wavelengths due to intrinsic spec- 
tral properties of the source and to the fact that at these wave- 
lengths flux contributions from an extended dust component 
complicates the determination of the flux density from Sgr A* 
(Eckart et al. 2006a Ghez et al. 2005b I. Concerning the X-Ray 
regime, there seems to be also one fairly firm case of QPOs in 
a flare observed by XMM (Belanger et al. l2006l l. 

A serious limitation for the possible detection of quasi- 
periodicities in NIR flares are observational constraints. 



Several factors must coincide: (i) the occurrence of a flare; 

(ii) the flare must be fairly bright and last long enough in or- 
der to be able to sample a sufficient number of oscillations; 

(iii) atmospheric conditions must be good (seeing ;S 0.8" and 
to ^ 3 ms at NACO/VLT) and - above all - stable during the 
entire observation, i.e. several hours. We estimate from our ob- 
servational experience that these conditions are fulfilled only 
10-20% of the total time dedicated to observations with NACO 
at the VLT. An additional problem is caused by gaps in the 
observational sequence that may be introduced into the time 
series due to the necessity of sky background measurements in 
NIR observations. Due to the extremely crowded field of the 
stellar cluster in the GC, the background cannot be extracted 
from the on-source observations. Also it is not clear whether 
each flare would be accompanied by a QPO. A new window 
into QPOs may have been opened by the discovery of Eckart 
et al. ( 2006b ), who found that possible periodic variations can 
be observed in polarized light while not being detectable in the 
overall flux. 

In the context of QPOs, it is also important to refer to the 
debate concerning the existence of QPOs in AGN. In many 
cases X-ray light curves of these objects appear to be charac- 
terised by red noise, which may lead to an overestimation of 
the significance of periodic signals in the light curves (see e.g. 
Vaughan 2003 Timmer & Koenig [T995l Benlloch et al. l2WlV 

In this section we re-analyse the so far best case for NIR 
quasi-periodicity in Sgr A*, that is the K-band flare from 16 
June 2003 which was originally reported by Genzel et al. 
(2003). Contrary to the other two mentioned flares, it has the 
advantage that apparent QPOs could be followed over 7 cycles. 
The imaging data were observed with the AO system/infrared 
camera NACO at the ESO VLT. The detector integration time 
was 10 s. Several images were averaged before saving the data. 
The effective time resoltuion is about one image every 40 s. 
We reduced the data in a standard way (sky subtraction, bad 
pixel correction, flat fielding). Applying a sliding window, the 
median of 5 images at a time was constructed. From these im- 
ages, the PSF was extracted via point-source fitting with the 
program StarFinder (Diolaiti et al. 2000 1. Subsequently, the 
images were deconvolved with a Lucy-Richard deconvolution 
and beam restored with a Gaussian beam. The flux of Sgr A* 
and of stellar sources in the field was extracted via aperture 
photometry with a sufficiently large aperture to contain the en- 
tire flux of the objects. Photometry calibration was done rela- 
tive to several stars in the field with a known flux. The result- 
ing extinction corrected light curve (assuming Ak = 2.8 mag) 
is shown in the upper panel of Fig.^ The black curve shows 
the flux of Sgr A* and the upper, green curve the flux of SI, a 
star at ~ 0.4" distance from Sgr A*. Its lightcurve is assumed 
to be constant. The lower, red curve shows the light curve of 
Sgr A* after subtracting a constant (0-40 min) plus a 5th-order 
polynomial fit (40-130 min) to the overall flare. 

The lower panel of Fig.^shows a non-oversampled Lomb- 
Scargle periodogram of the Sgr A* flux after subtraction of the 
overall flare. There is a broad peak visible around 2 x 10~ 3 Hz 
that corresponds to the 16-21 min variability that can be seen 
in the light curve. In order to examine the significance of this 
peak, we followed the recipe of Vaughan (2005 1. The assumed 
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Fig. 1. Sgr A* NIR flare from 16 June 2003. Upper panel: Light 
curve of Sgr A* in black. A polynome fit to the overall flare and 
the residuals after subtracting this fit are shown in red. The light 
curve of a constant star, SI, is shown in green. Lower panel: 
Periodogram of the Sgr A* flare after subtraction of the over- 
all flare, i.e. corresponding to the lower, red curve in the upper 
panel. The straight red line shows a de-biased power law fit to 
the data points when considering all frequencies greater than 
2 x 10~ 3 Hz. The dashed lines give the 1-, 2-, 3- and 5<x confi- 
dence bands derived with the method of Vaughan ( 32005> . 



underlying power law is indicated by the straight red line in 
the lower panel of Fig. ^ This power law with an index of 
-0.99, i.e. a classical red noise, was found by fitting all data 
points in the frequency range > 2 x 10~ 3 Hz. The prescription 
of Vaughan ( 2005 1 takes the bias into account that originates 
when fitting the logarithm of the periodogram. Dashed lines 
in the plot of the periodogram indicate the 1-, 2-, 3- and 5<x 
confidence bands which can be calculated analytically. We also 
used Monte Carlo simulations to derive these confidence limits 
which leads to the same result. 

It can be seen that the peak around 10~ 3 Hz reaches the 5<x 
threshold. Even after correcting for the model uncertainty, as it 
is suggested in Vaughan ( 2005 1, the significance is still at the 
4.2cr level. 

The critical point in this procedure is the choice of the 
slope of the power spectrum. We followed the prescription of 
Vaughan ( 2005 1. We think our approach is justified because the 
periodogram shows a clear component of red noise (see also 
discussion in Belanger et al. 2006 1. 



We conclude that, on the one hand, there are strong indi- 
cations for a quasi-periodicity in the NIR flare from 16 June 
2003. On the other hand, this single flare by itself is not enough 
to present ironclad evidence. However, two more events were 
reported that support the presence of a periodicity in the con- 
sidered time range: the second K-band flare reported by Genzel 
et al. (2003 1 and the flare observed in polarimetry by Eckart et 
al. J2006b> . In fact, everytime sub-flares have been observed so 
far, they showed a period that lies within the 20+3 min interval. 
The likelihood of seeing sub-flares at this interval in consecu- 
tive observations is small. The probabilities multiply and can 
be estimated simply: the probability to detect a peak in the pe- 
riodogram of a single observation is at most at the percent level, 
see discussion above. That means that the chance to observe a 
peak at the same frequency interval in three consecutive runs is 
~ 0.1%, and ~ 0.01% if a fourth NIR flare showing sub-flares 
has a period at ~ 20 min. Therefore we believe that the evi- 
dence for an intrinsic cause of these modulations must be taken 
serious. 

3. The Hot Spot/Ring Model 

Disk accretion is generally thought to be the source of galac- 
tic nuclear luminosity and the driver of activity whenever the 
angular momentum dominated disk-type accretion takes place. 
In the case of Sgr A*, the accretion rate onto the central black 
hole and the net angular momentum of the accreted gas appear 
to be very small 10~ 7 Mq yr~\ see Bower et al. 2003 ) and a 
persistent accretion disk is not present (Coker & Melia 1997 1. 
Despite this fact, the very occurrence of the flares and the 
minute time-scale of their rise and decay indicate that the im- 
mediate neighbourhood of the black hole is not entirely empty; 
likely, episodic accretion takes place. In fact, Moscibrodzka et 
al. |2006 suggest that a transient torus develops at a distance 
of a few or a few ten gravitational radii and exists for several 
dynamical periods. The inner disk/torus is subject to violent in- 
stabilities, but it is difficult to identify the dominant mechanism 
responsible causing its destruction. This makes it tempting to 
attribute the outbursts to some type of disk instability, for in- 
stance of the limit-cycle type as known in dwarf novae (albeit 
with the possibility of a physical cause different from the hy- 
drogen ionization/recombination there). Disks in the vicinity of 
black holes live in a relativistic environment, and proper mod- 
eling requires relativistic (magneto-)hydrodynamics (see e.g. 
Tagger & Melia 2006 1. In the following, however, disk insta- 
bilities can be ruled out by such a wide margin that a non- 
relativistic order-of-magnitude treatment suffices our purpose. 

So close to the black hole, self-gravity will not play a role, 
allowing for the standard one-zone approximation of the verti- 
cal hydrostatic equilibrium, which relates the azimuthal veloc- 
ity v^, the sound velocity c s , the (half-)thickness of the disk h, 
and the radius s by 



The rise time Af + of a disk outburst may be estimated by 
the radial range of the disk involved in the outburst, As, and the 
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propagation speed of the outburst, which is ac s (Meyer 1984 
see also Liu & Melia l2"0"0"2"l ). 

As As As s s As 

Af + ~ > — = — - ~ fdyn 7 — (3) 

a c s c s h h s 

with the dynamical timescale fd vn ~ — ■ Here a is a parameter 
related to the disk viscosity and is of order unity or smaller. For 
the outburst to be a disk phenomenon and not just something 
happening in a small ring, it is required As » s. 

Even for an only moderately thin disk of, let's say, h/s = 
1CT° 5 this then requires rise times which are more than an order 
of magnitude longer than the dynamical time scale. The small- 
est dynamical time scale allowed for is the one at the ISCO. As 
discussed above, however, the entire outburst (rise plus decline) 
encompasses only a few dynamical time scales. 

This ratio of the observed time scales makes a normal disk 
instability an exceedingly unlikely cause of the outbursts. It fur- 
thermore shows that, at best, a relatively small radial range of 
the disk {As ~ s), close to the ISCO, may be involved in what- 
ever causes the outbursts. 

It is, however, noteworthy that, for a moderately thin ac- 
cretion disk, the duration of the outburst is of the order of the 
timescale a sound wave needs to travel around the black hole 
once. This makes a scenario attractive in which the outburst is 
caused by some non-axisymmetric process (be it a disk insta- 
bility of a non-axisymmetric type, confined to the innermost 
disk regions, be it material falling into the disk there locally), 
and a sound wave originating from this event and traveling the 
inner disk until the outburst comes to an end when the actual 
source of it has been accreted into the black hole. While the 
sub-flares tracing the actual local event disappear suddenly as 
soon as that region has been accreted, the outburst itself ceases 
somewhat slower as the soundwave, of course, travels not only 
azimuthally, but also radially outwards, though only a small 
distance as shown above. 

Motivated by these considerations, we adopt a two com- 
ponent hot spot/ring model in order to fit the data presented 
in Eckart et al. ( 2006b I. In this model the broad underly- 
ing flare is caused by a truncated disk/ring. The sub-flares 
are due to a compact emission region on a relativistic orbit 
around the MBH. Such hot spots can be created in reconnec- 
tion events, in which thermal electrons are accelerated into a 
broken power law (Yuan et al. |2004). General relativistic ef- 
fects imprint on the synchrotron radiation of such inhomo- 
geneities (Broderick & Loeb l2005ll2006l Dovciak et al. 120041 
Hollywood et al. [T995l Hollywood & Meha flW71 Bromley et 
al. l2U0Tl Falcke et al. l2000l Pineault [T98"TV Redshifts, lensing, 
time delay, change in emission angle and change in polarization 
angle belong to these. The code by Dovciak et al. ( 2004 1 takes 
these special and general relativistic effects into account by us- 
ing the concept of a transfer function (Cunningham 1 1 9751 1 . A 
transfer function relates the flux as seen by a local observer co- 
moving with an accretion disk to the flux as seen by an observer 
at infinity. This transfer of photons is numerically computed by 
integration of the geodesic equation. For the change of polar- 
ization angle the method of Connors & Stark ( 1977) has been 
used. In geometrical optics approximation, photons follow null 



geodesies and their propagation is not affected by the spin-spin 
interaction with a rotating black hole (Mashoon 11973b . This 
means that wave fronts do not depend on the photon polariza- 
tion, and so the ray-tracing through the curved spacetime of the 
black hole is adequate to determine the observed signal. 

Following the above discussion we have fixed the radial ex- 
tent of the disk to 2R$ , beginning at an inner edge. We realize 
that there is a range of definitions for an inner edge that depend 
on the accretion rate and may be non-axisymmetric and time- 
variable (Krolik & Hawley 2002). They are, however, generally 
located in the vicinity of the marginally stable orbit. For sim- 
plicity we have assumed that for a* S 0.5 both coinside, and 
our general result - that the observed time variability is com- 
patible with a discription of a spot/ring combination orbiting a 
central MBH - does not depend strongly on this assumption. 
Due to the observed timescale of the QPOs, the spot is within 
the marginally stable orbit for a* S= 0.5. In this case we as- 
sumed the spot to be freely falling and the disk to have its inner 
edge at ~ Rs /2 above the event horizon. Magnetic field effects 
could bar the spot from freely falling (Krolik & Hawley 2002 1 
but require full relativistic MHD simulations that are beyond 
the scope of this work. 

The disk/ring's time behaviour was assumed to be Gaussian 
and to account for the main flare. The spot is orbiting within 
the disk and the equatorial plane of a Kerr BH. Its intrinsic lu- 
minosity follows a two-dimensional normal distribution with 
cr ~ 1 .5Rs , but being cut off in the radial direction to fit within 
the disk. The spot follows a circular trajectory (a* ^ 0.5) cho- 
sen such that the observed periodicity is matched. Note that the 
radius for the corresponding orbit is spin dependent, see equa- 
tion or Bardeen et al. (119721 . In the case a* £ 0.5 the spot 
is freely falling and it has to be checked whether the timescale 
can be matched within its uncertainty. 

We identify all three sub-flares with the same spot. This can 
well be matched with the corresponding synchrotron cooling 
timescale (Eckart et al. l2006al b: Yuan et al. EOOSlEOOfl . The 
observed differences in the flux of each sub-flare are very likely 
due to intrinsic changes of the spot, i.e. the spot is very likely to 
evolve. Regarding the fact that the physics of such an inhomo- 
geneity in the accretion flow is unknown at present, we simply 
model the changes of the intrinsic luminosity by hand. This is 
done by changing the value for the spot luminosity discontinu- 
ously instead of some unphysical parametrization. Normalized 
to the third revolution (~ 36 - 54 min) the spot has an intrinsic 
luminosity of 50% between 0-18 min, 20% between 18-36 min 
and 10% after 54 min. The small jumps seen in some light 
curves are an artifact of this procedure. 

Another approach would be to assign different spots to dif- 
ferent sub-flares. Assuming a typical spot-lifetime of approx- 
imately one revolution and that every spot is on roughly the 
same radius, multiple single spots can reproduce a peak in the 
periodogram (Schnittman 2005). This procedure could natu- 
rally lead to better fits, because the phase of each spot would be 
a free parameter for every single sub-flare. We rejected that ap- 
proach, however, as it only increases the already large number 
of degrees of freedom. 

To investigate the parameters of Sgr A* we have fitted for 
the inclination angle, the dimensionless spin parameter a*, the 
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Fig. 2. Sketch of the Newtonian behaviour of the ^-vector (red) 
for the two magnetic field configurations. The black circle rep- 
resents the MBH, the ellipse within the equatorial plane indi- 
cates a projected Keplerian orbit of a hot spot in the local co- 
moving frame. Left: the case where the £-vector is always con- 
stant. Right: the case of an azimuthal magnetic field (like the 
ellipse) where the £-vector of the emitted synchrotron radia- 
tion rotates. 

brightness excess of the spot with respect to the disk, the initial 
phase of the spot on the orbit, the orientation of the equatorial 
plane on the sky and the polarization degree of the disk (con- 
strained to < 15%) and the spot (< 70%, a value that can well 
be achieved via synchrotron radiation). As the last five param- 
eters can be changed after the ray tracing, it costs little com- 
putational effort to find the least x 2 values for given spin and 
inclination (Lourakis 2004 1. For these two parameters we have 
chosen a discrete grid with < a* < 1 and steps of 0.1 and 
an inclination angle 10° < i < 70° with steps of 5°; i = 0°: 
face-on. For i > 70° multiple images could become important, 
which are not included in our treatment (see next Sect.). 

As the observational data include the effects of foreground 
polarization, we also take these into account. Following the cal- 
ibration of Eckart et al. J2006b> we fix it to be 3.4% at 25°. 
As only rough estimates can be given on the magnetic field 
of the spot, we did the fits for two different field configura- 
tions as an approximation (Pineault 19^}. The first configura- 
tion is such that the resulting projected E- Vector is always per- 
pendicular to the equatorial plane. Such preferred orientation 
could result from perturbations in the disk similar to sunspots 
(see also Shakura & Sunvaev fl973> . As a second configuration 
we have allowed for a global azimuthal magnetic field. This 
behaviour may be caused by the magneto-rotational instabil- 
ity and is motivated by global MHD simulations (Hirose et al. 
2004 1. Contrary to the first case, this configuration leads to a 
rotation of the ^-vector along the orbit from a Newtonian per- 
spective (see Figure 2). 

4. Results & Discussion 

Figure[3]shows the fit with the least x 1 value within the discrete 
grid discussed in Sect. [3] in case of a constant ^-vector. High 
inclination and high spin give the best solutions. The case of 
the azimuthal magnetic field can be seen in Figure|4] Here high 
inclination and medium spin is preferred. Within the hot spot 
model the degree of polarization does not necessarily have to 
be fitted. The spot may have influence on the disk and may 
polarize the disk flux. 

The confidence contours shown in Figure |5]reveal that for 
both magnetic field cases the spin paramter a* cannot be well 
constrained from the current data with its uncertainty. Within 
the 3<x limit a* can be inferred to lie in the range 0.4 < a* < 1. 
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Fig. 3. The best fit solution (in red) for a constant E- vector. 
Shown is the flux (top), polarization angle (middle), and po- 
larization degree (bottom). The parameters of the model are 
i = 70°, a* « 1, a location of the projection of the disk on the 
sky of 100° and a polarization degree of the disk (spot) of 4% 
(60%). The spot is orbiting at a radius r - 4-.4r g . 

The range of the inclination is different for the two magnetic 
field scenarii. As the ^-minimum is lower for the constant E- 
vector, this case should be weighted more. That means the in- 
clination is 35° on a 3<x level. 

Our solution with a* > 0.4 is in agreement with Genzel et 
al. (2003 1 and Belanger et al. (2006), who inferred values of 
a* > 0.5 and a* > 0.22, respectively. They arrived at these 
lower limits by interpreting the peak in the power spectra of 
their measured light curves using the hot spot model. As it is 
not possible to decide where the spot is orbiting, only a lower 
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Fig. 4. The best fit solution (in red) for a global azimuthal mag- 
netic field. Shown is the flux (top), polarization angle (middle), 
and polarization degree (bottom). The parameters of the model 
are inclination = 70°, a* ~ 0.5, a location of the projection of 
the disk on the sky of 0° and a polarization degree of the disk 
(spot) of 9% (60%). The spot is orbiting at a radius r = 4.5r g . 

bound can be given with the equality applying if the spot is 
exactly on the ISCO. Only polarimetric observations together 
with a relativistic modeling can improve on this within the spot 
model due to the additional indepent information. Our analysis 
of recent polarimetric measurements, however, shows that very 
accurate data are needed to give strong constraints on the spin 
parameter. 

One should be careful about the role of general relativis- 
tic effects for polarization fluctuations and there does not seem 
to be a clear way towards a unique interpretation at the present 
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Fig. 5. The confidence contours for the constant ^-vector (top) 
and the azimuthal magnetic field case (bottom). The red (green) 
line is chosen such that the projection onto one of the parame- 
ter axes gives the lcr (3cr) limit for this parameter. The respec- 
tive ^-minimum is marked by the big cross. Note the different 
scales on the abscissae. 



stage. Likewise the inferred values of a* and the inclination are 
subject to many uncertainties - mainly concerning the geome- 
try of the source and its intrinsic polarization. 

We note that we did not take higher order images into ac- 
count. While they offer the unique possibility of testing general 
relativity in a strong field quantitatively (Broderick & Loeb 
2006 1, their detection is extremely difficult and lies well be- 
yond present observational capabilities. Interferometry in the 
NIR or sub-mm domain may resolve this question in the fu- 
ture. Only for special geometries, where strong lensing is very 
important, they yield a non-negligible contribution (Broderick 
& Loebl2006l Bozza & Mancini 2004J Horak & Karas l2006al 
2006b I. The lack of multiple images in our models limits their 
applicability to i % 70°. One may, however, speculate, what 
signatures of multiple images to look for. Time delay between 
1st and 2nd images is a very characteristic number, related to 
BH mass: it is given by the circumference of the photon circu- 
lar orbit, R p = 2(1 + a)GM/c 2 , where a = cos[| arccos a*] 
(Bardeen et al. 1972). Hence, the expected delay is of the or- 
der of 575 (M/(3.6 x 10 6 M o ) sec in the Schwarzschild case 
(a* = 0; this interval becomes roughly a factor of 3 shorter 
in case of a maximally rotating black hole and prograde ori- 
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entation of the photon trajectory, while it is about 1.3 times 
longer for the retrograde orientation). Unfortunately, the sec- 
ond image provides an order of magnitude less photons, and 
higher order (indirect) images are still weaker because the ob- 
served radiation flux decreases exponentially with the image 
order n increasing (its observed flux is roughly proportional to 
exp(-27rn); cf. Luminet 1979 1, unless, again, a special geome- 
try is assumed in which caustic lensing occurs. 

In fact, the importance of significant lensing enhancement 
of the observed flux from accretion disks has been tradition- 
ally neglected, because it would require large inclination an- 
gles. Such geometry is not likely in case of standard (optically 
thick) accretion disks, because self-obscuration does not allow 
to observe the indirect images. However, the situation is differ- 
ent for optically thin flows (as noticed by Bursa et al. 2004 1. 
Therefore, Sgr A* is a suitable object in which the flux fluctua- 
tions could be considerably modulated by strong-gravity lens- 
ing if the inclination of the putative disk is very large (see also 
Falcke et al. 120001 Bromley et al. l2"00TV 

In case of non-zero rotation of the BH, the expected 
lightcurves exhibit stronger fluctuations because the struc- 
ture of light-ray caustics is more complex (e.g. Hollywood & 
Melia HW71 Rauch & Blandford 1994 Viergutz|T333}. Rapid 
decay of the signal strength with increasing n holds also for 
a rotating black hole, although the prospect of detecting the 
higher-order images appears even more tempting, because their 
mutual delays at the point of arrival and relative polarization 
degree can set tight constraints on the black-hole angular mo- 
mentum. In conclusion, higher order images are not possible to 
measure yet, though it would be very useful and should be at- 
tempted in the future because their interpretation could be less 
ambiguous. 

5. Summary 

After showing that the quasi-periodicity of the June 2003 flare 
(Genzel et al. 2003 1 appears to be clearly significant, we have 
outlined a consistent physical picture for the NIR flares of 
Sgr A*. In that model the broad flares with a typical timescale 
of 60 - 100 min are due to a sound wave that travels around 
the MBH once, caused by a non-axisymmetric perturbation in 
the innermost region of an accretion disk. The sub-flares su- 
perimposed on the broad flare are due to an orbiting compact 
source of transiently heated and accelerated electrons (Yuan et 
al. |2)04 1, whereas the acceleration mechanism may be linked 
to the local event causing the sound wave. As we have shown, 
modeling this with the simple hot spot idea leads to reasonable 
fits of the polarimetric NIR observations presented by Eckart 
et al. (2006b I. The least x 1 values are obtained for a spin pa- 
rameter a+ close to unity, on a 3<x level it is constrained to the 
range 0.4 < a* < 1. The inclination tends to be high, it can be 
constrained to i > 35°. 
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